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We discuss the dynamics of a Bose-Einstein condensate of atoms which is suddenly coupled to a 
condensate of molecules by an optical or magnetic Feshbach resonance. Three limiting regimes are 
found and can be understood from the transient dynamics occuring for each pair of atoms. This 
transient dynamics can be summarised into a time-dependent shift and broadening of the molecular 
state. A simple Gross-PitaevskiT picture including this shift and broadening is proposed to describe 
^^ the system in the three regimes. Finally, we suggest how to explore these regimes experimentally. 
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H I. INTRODUCTION 

^H Optical and magnetic Feshbach resonances [H [2] occur when two colliding particles are resonant with a bound state 

V/^ [3]. They have been extensively used to probe and control the properties of ultracold atomic gases [4]. In particular, 

they give the ability to tune the interaction between atoms, making ultracold gases very useful to test many-body 

"TT theories for a wide range of interaction strengths, which is generally not possible in other kinds of systems. These 

D resonances also led to the controlled association of atoms into diatomic molecules, a process called phot oassociat ion 

F^ [5J in the case of an optical Feshbach resonance. This association is a coherent process and when applied to degenerate 

Q gases, it can result in the formation of a Bose-Einstein condensate of molecules, also called molecular condensate [6l[7]. 

j ! In this work, we discuss the dynamics of an atomic Bose-Einstein condensate which is suddenly coupled to a 

Cd molecular condensate through an optical or magnetic Feshbach resonance. In the case of an optical resonance, this can 

Q be achieved by turning on a laser on resonance between a colliding pair of atoms and an excited bound state. In the 

I case of a magnetic resonance, this can be achieved by suddenly imposing a magnetic field which brings a bound state 

^^ in resonance with the colliding pair. In both cases, the two-atom system can be effectively described by two channels: 

b^ the (open) scattering channel where atoms collide through an interaction potential U{r) and a (closed) molecular 

Q channel where the atoms can bind through an interaction potential Um{^)- The vector r denotes the relative position 

•—I of the two atoms. When r becomes large, the potential U{f) -^ , and the potential Umir) goes to the internal energy 

Em- The two channels are coupled by some potential W{f). In the case of an optical resonance, the two channels 

correspond to different electronic states and the coupling is due to the laser. In the case of a magnetic resonance, 

1^ the two channels correspond to different hyperfine states and the coupling is due to hyperfine interaction. With this 

^^^ simple picture of the resonance at the two-atom level, we can write equations for a Bose-Einstein condensate. 

m 

o 

] II. EQUATIONS FOR COUPLED ATOMIC-MOLECULAR CONDENSATE 

r^^ A coupled atomic-molecular Bose-Einstein condensate can be described by time-dependent equations for its cumu- 

^^ lants. Truncating to first order as prescribed in Ref. [8j, one obtains the following set of equations [9l[T0l[TI]: 

^ ih—^\x,y) = {H^^Hy)^\x,y)^U{x-y)^{x,y)^W{x-y)^m{x,y) 

d 

ih—^m{x,y) = {H:,^Hy-ij/2)^rn{x,y)^Um{x-y)^m{x,y)^W{x-y)^{x,y) 

where x and y are space coordinates, t is the time variable (for clarity, it is omitted but implicitly present as an 

argument of all functions), h is the reduced Planck's constant, H^ = — ^jf- + V{x) is the one-body Hamiltonian for 
a single atom (M is the atom mass and V some external trapping potential), ^ is the condensate wave function, ^^ 
is the pair wave function in the closed channel, and $ is the pair wave function in the open channel. The function 
^'{x^y) = ^{x^y) — ^{x)^{y) is the pair cumulant. It corresponds to the deviation of the pair wave function from 
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its noninteracting form, which is a product of two condensate wave functions. We included a loss term —i^/2 in the 
closed channel to account for possible decay by spontaneous emission. 

Usually, only one molecular level (frn of the potential Um is in resonance with the atomic condensate pairs. Thus, 
we can write: 

^^{x,y,t) = ^m{R.t)^m{r) (2) 

where R = {x-\-y)/2 and r = ^— x are the centre-of-mass and relative coordinates. ^^ is a one-body field corresponding 
to the molecular condensate wave function, describing the motion of a molecule whose internal state (pm satisfies 

( - ^ + Um{r))^m{r) = {Em - E^)^m{T) 

where £^5 is the binding energy associated with the molecular state iprn- We choose the normalisation {(fm\^m) = 1, 
and call A = Em — £^5 the detuning from the resonance. 

After some approximations detailed in the appendix, we find that ^ and ^m satisfy the closed set of equations: 



ih^{x) = (-^ + n^)+^|^(f)p)^(x) + ^*(x)(l+^a/(t))^^^(f) (3) 



ih 



5*„^_. / tl'V 



dt 



(f) = (-^+2na*)+A-A'-^7/2 + w;V(i))*n^(a*) + w'*(a*)' (4) 



where g is the coupling constant for elastic collisions between atoms, w is the coupling constant for atom-molecule 
conversion, and A' is a light shift of the molecular level due to the coupling. For three-dimensional systems, we have 

a = 7r/2 

The term w'^afit) can be interpreted as a time-dependent shift and broadening of the molecular level. It is a time- 
dependent version of the static light shift A^ when the interaction is switched on, a large shift and broadening of the 
molecular level appears during some transient regime and then eventually goes to the static value A^ The reason for 
the broadening is the loss of molecules into atom pairs. Indeed, during a certain time, the coupling w can dissociate 
the molecules into the pair continuum instead of dissociating them back to the atomic condensate state, a process 
called "rogue dissociation" in Ref. ^12j . 

III. COHERENT, ADIABATIC AND DISSOCIATION REGIMES 



We now investigate the different regimes that one can get from Eqs. (|3j|4|. These regimes can be identified by the 
qualitatively different short-time dynamics. At short times, we can neglect the motion of molecules in the trap. We 
can also write ^^^(x, t) ^ {^rn{x^t) — ^^(x, 0))/t. Using this with the initial conditions (26) in Eq. (H), we can 



solve for ^m(^) and insert it into Eq. (|3|. We obtain an effective Gross-Pitaevskii equation for the condensate wave 
function ^ 

^^?(-i = f - ^ + nx) + (9- w^ '-^^ ^) i*(f)|2)*(f) (5) 

dt^ ' \ 2M ^ > \y A - A' - i7/2 + u'2a/(i) - f / ' ^ >^ J ^ ^ ^^ 

In the local density approximation, this equation leads to the familiar rate equation 

p{x,t) = -K{t)p{x,tf (6) 

where p{x,t) = |*(a;, t)p is the atomic condensate density, and K{t) is a time-dependent rate coefRcient 

^ ^ h A - A' - i-f/2 + w^af{t) -f ^ ^ 

This rate coefficient depends only on molecular physics parameters g, w, A', and 7. In fact, it can be shown to 
derive from the time-dependent two-body theory, except for the term proportional to f{t) in the numerator of K{t). 
We find that this term can be neglected when we compare with numerical calculations based on Eqs. ( 19|2Q|25 ). This 



indicates that the term f{t) should be discarded from Eq. (|3|. Depending on the relative strengths of the terms in 
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Figure 1: Schematic evolution of the instantaneous rate coefficient (thick black curve) showing the three regimes (a), (b) and 
(c). Upper graph: mode rate co upling (i^ ^ t^)- Lower graph: strong coupling (iw ^ Ia)- The grey curves correspond to 
the limiting expressions ( 8|9|1Q ). Note that the parabolic curve associated with Eq. (lo]) is a universal upper limit of the rate 
coefficient. 



the denominator of K(t\ the rate coefficient for a fixed detuning goes through three subsequent regimes illustrated in 
Fig. IT] linear with time (a), square root of time (b) and constant with time (c). 



4. / h\^/^ 



Kit) = -[j^j Vi ior t^^KtA 



K{t) = -Im 



where we define the two-body time scales 



h A - A' - ij/2 



hV n^' 



for tA < t 



(8) 

(9) 

(10) 



(7r/2)2 \MJ \w 
tA = (V2)'z^hr 



W^ 



h?\h J (A - /^'Y + (7/2)2 

The first two regimes correspond to transient dynamics of the pairs, while the third regime yields the rate coefficient 
of the time-independent two-body theory. Note that the second regime might not exist for small coupling w^ when 
't'w > tA- Then the first and third regime are separated by the time scale tj = ^/tw^A- 

In principle, the two transients regimes (a) and (b) always occur, provided the coupling is large enough. However, 
they might often occur at very small times, so that they do not affect the system significantly. In other words, the 
time scales tyj and Ia might be much smaller than the time scale for the evolution of the gas. We define this time 
scale as the time tp needed to reduce the atomic condensate population by a half. Obviously, this time scale depends 
on the initial density po of the gas, as can be seen from Eq. (6|. Three situations are possible: 



• tp ;» t^, t^, in which case the evolution of the gas is governed by the constant rate coefficient (10) 

• t^ <C tp <C t^, in which case the evolution of the gas is governed by the transient rate coefficient ^ 

• tp <^tyj^t^^\VL which case the evolution of the gas is governed by the transient rate coefficient (|8| 
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Figure 2: Different regimes of Feshbach resonance in a condensate as a function of the detuning and width ^^(A — A^)^ + (7/2)^ 
and the coupUng strength w^ , for a fixed density. The boundaries are indicated by dashed Unes and equaUties of times scales 
defined in the text. 



We call these three situations the adiabatic regime, the dissociation regime and the coherent regime, respectively. 
They are re presented in Fig, [g] for a fixed initial density, as a function of the coupling strength w'^ and "width and 
detuning" y^(A — A')^ + (7/2P. In the adiabatic regime, the molecular condensate ^^ has a very small population 
because it decays too fast or is too far from resonance and it can be eliminated adiabatically. The system is then 
described at all times by a single Gross-Pitaevskii equation such as ([5|. In the dissociation regime, the molecular 
condensate has also a small population, but now this is because it is dissociated by the coupling into the pair continuum 
[T2] . The system can be described only in terms of the atomic condensate and dissociated pairs. Interestingly, it is a 
universal regime in the sense that all quantities depend only upon the mass of the atoms, and not on the particular 
resonance used. In the coherent regime, the molecular condensate can have a large population and dissociates only 
into condensate atoms. This creates coherent oscillations between the atomic and molecular condensates, which are 
described at all times by two coupled Gross-Pitaevskii equations. This regime corresponds to the superchemistry 
described in [13j. 

We stress that in all cases, the short-time dynamics of the atomic condensate is described by the rate equation 
(|6| with the time-dependent rate coefficient Q . This shows that the loss of condensate atoms at short times goes 
like t, t^/^ and t^ for respectively the adiabatic, dissociation, and coherent regimes. For longer times, the system is 
well described by the two Gross-Pitaevskii equations (|3][4| containing the transient shift and broadening. However, 
higher-order cumulants (quantum fluctuations) neglected in ([l]), as well as inelastic collisions fT4] between molecules, 
dissociated pairs or condensate atoms play a role for long times (typically ~10 /iS in the cases shown in Fig. [sj. 



IV. APPLICATIONS 



Experimentally, the adiabatic regime has been well explored with both optical and magnetic Feshbach resonances 
dSl [161 [T71 [E] . As explained above, the system is then described by a Gross-Pitaevskii equation where the scattering 
length appearing in the mean-field term is changed by the resonance. In the case of an optical resonance, the scattering 
length has an imaginary part because of losses from the molecular state by spontaneous emission. These losses are 
described by the rate equation ([6| with the constant rate coefficient (10). 



The dissociation regime has proved more difficult to observe. For optical resonances, the losses from the molecular 
state usually confines the system to the adiabatic regime, as can been seen from Fig.jS] To reach the dissociation regime, 
one needs to increase the coupling, ie the intensity of the laser. A high-intensity experiment was performed at NIST 
with a sodium condensate ^6j, but the laser power was not sufficient to reach the dissociation regime. Experiments 
in the dissociation regime were performed with magnetic resonances by switching the magnetic field suddenly on or 






Figure 3: On-resonance photoassociation (A — A^ = 0) of a uniform atomic condensate of sodium (a), ytterbium (b) and 
strontium (c), for typical transitions. In each row, the right panel the time evolution of the populations, according to Eq. 
( |19|20|24' ), based on the paramaters indicated by the black dot in the left panel, which is a regime diagram similar to Fig. [5] 
Solid curve: atomic condensate population; dotted curve: dissociated pair population; dot -das hed curve: molecular population. 
The short-dashed curve shows the atomic and molecular population following from Eqs. (3|4). 



We used the following parameters: w = 1.255 • 10"^^ J-m^/^ -f/h = 18 MHz, w = 3.850 • 10"^^ J-m^/^ -f/h = 364 kHz for case 



(b), 



: 6.347 • 10 ^° J-m^/^, ^y/h = 15 kHz for case (c). For all three cases, the initial density is po = 6 • lO^^cm 



very close to resonance [19j. They resulted in an explosion of hot atoms, and thus were called "Bosenova" experiments. 
The hot atoms have been identified by several authors with the dissociated pairs discussed above [20l[2T], although 
several alternative theories exist [22j |23l |24l [25] . Quantitative comparison with any theory has not been completely 
satisfactory so far [26j, suggesting that either some element is missing from the theory or some experimental condition 
has been misunderstood. 

To our knowledge, the coherent regime has not been observed. Here we propose to observe the dissociation and 
coherent regimes by using optical resonances in systems with narrow intercombination lines, for which the molecular 
states are long-lived. Figure shows the cases of ytterbium and strontium, for typical resonances. In the case of 
ytterbium, it appears possible to reach the universal regime of dissociation for laser intensities of about 3 W/cm^. 
In the case of strontium for similar intensities, one period of coherent oscillation between the atomic and molecular 
condensates can be observed. 

To illus trate some approximations made in the previous section, we made numerical calculations based on equations 
( 19|20|24 ) and (3J4) in the case of a uniform system. Fig. [s] shows the evolution of the populatio ns in the atomic, 
molecular condensate and dissociated pairs. One can see that the short-time dynamics of equations J^ 19|20|24 ) is well 
described by the coupled Gross-Pitaevskii equations with the transient shift and broadening, Eqs. (|3][4|. 

Finally, we present in Fig.|4]the energy distribution of the dissociated pairs for the different systems. In the universal 
regime, the dissociated pair distribution can be calculated analytically 
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(11) 



where Erf is the error function. From this, one finds the energy distribution p{E]^) ex -^2k\Ck\'^. Figure H shows the 
average energy of the dissociated pairs as a function of time, and the final energy distribution. Althougn~Eq. (11) is 



expected to work only at short time and in the universal regime, it gives a good qualitative prediction of the energy 
distribution (up to a normalisation factor) . It is worth noting that when the laser is turned off quickly, although the 
number of dissociated pairs is maintained, their energy might go down. We find that the laser has to be turned off 
very fast to maintain the final energy distribution. Therefore there is a range of switch-off speed for which the number 
of pairs is maintained, while their energy is lowered. 



V. CONCLUSION 



We showed that three limiting regimes are expected when coupling an atomic condensate to a molecular condensate 
by an optical or magnetic Feshbach resonance. The three regimes can be identified according to the short-time dynamics 
of the system. This short-time dynamics can be understood at the two-body level by the transient response of the 
atom pairs, which is well described by a time-dependent shift and broadening of the molecular state. This shift and 
broadening can be incorporated into a set of coupled Gross-Pitaevskii equations to describe the short-time dynamics 
of the system. Finally, we suggested different systems to observe the three regimes experimentally, in particular the 
universal regime where atom pairs dissociate, and predicted their short-time energy distribution. 
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Figure 4: Energy of the noncondensate atoms (dissociated from the molecular state and not returning to the condensate) for 
the cases (b) and (c) of Fig. Is] (ytterbium and strontium). Note that all graphs are in logarithmic scales. Left panels: average 
energy as a function of time. For long times, the average energy decreases as 1/t curves - the dashed lines are arbitrary oc 1/t 
curves shown for reference. Right panel: energy distribution at t = 5/iS for case (b) and t = 10 /iS for case (c). The dashed 
curves are obtained from Eq. (11) . The jagged appearance of the curves comes from the numerical discretisation of the pair 
continuum. 



VI. APPENDIX 

In this appendix, we give the derivation of Eqs. ( 3p ). 



Introduction of coupling constants 

We first express the pair cumulant <!>' in terms of centre of mass and relative motion by expanding it over the 
eigenstates ^^(r) of the relative Hamiltonian 



2v72 



M 



/7(r))^,Kr) = ^,^^,Kr), 



where Er = h^k'^ /M. One can write: 



$'(f,y,t) = ^'{R,?,t) = j ^^Cg(i?,t)(pg(r). 



(12) 



For simplicity, only the continuum eigenstates </5g(r) of the spectrum are taken into account (the bound states are 

supposed to be far off-resonant) and we choose the following normalisation: (^'gl'/Jp-) = {2t:)^ 5{k — p) . 
Inserting the decompositions (J2| and (12) into the equations (jTl), one obtains: 



a* 



d_ 



— (f) = iJ,*(f) + g{x) + j **(f + x) {U{r) j j^C^ui\ + ^)ni^ + W^(0*™(^ + x)vUr))d?r 



if^ilC^uiR) = i-^+Eu)Cm 



4M 






where 



g{x) = / **(f+f)C/(r)*(f + -)*(f - -)d^f 

gj:{R) = j ^^^{r)U{r)^{R+'-)^{R-'-)<fr 

w{R) = j ^m{^W{r)-^{R+'^)-^{R-'^)cPr 

Vm{R) = Jvm{^[v{R+^) + V{R-'^)'^^^{r)d^f 

%-(^) = / n(^ {viR + ^) + V{R - ^))<p,-(rlrf 



¥ 



are various matrix elements of the couplings U^W, and V between the states ^, ^^ and C^. 

The condensate wave function ^ and centre-of-mass function ^^ and C^ extend over mesoscopic scales, typically 
up to the micrometre range. On the other hand, quantities such as the interaction potential [/, Um and the molecular 
state (pm extend over microscopic scales, on the nanometre scale. Therefore, we can make the following approximations 



g{x) ^ ( f U{r)d^A^%x)^{x)^ =g^%x)^{x) 



Vm{x) ^ 2V{x) 



which leads to 



ih—{x) = if,*(f) + vl/*(f) (g^{xf + j ^^5-q-(f) + w^<i>Ux)) (13) 



dt ' ' V 4M ' ' " J "'' ' ' ' J (27r) 

This is not the most convenient representation for this set of equations, because in the case of a potential U{f) which 
is strongly repulsive at short distance, quantities such as g might diverge. Often, one treats this problem by replacing 
the actual potential U{f) by a contact pseudopotential g5^{f), where g is adjusted to some relevant finite value. This 
formally replaces g hj g, but also introduces ultraviolent divergences, due to the infinitely short range nature of the 
delta function. One then has to renormalise the coupling constant g using standard techniques in field theory, so as to 
reproduce physical results. However, in the present case, the whole set of equation is in fact well-behaved with respect 
to the interaction U and we do not need to resort to such renormalisation techniques, but simply rewrite the equation 
in another representation. 

We decompose C^ into two contributions, the adiabatic and dynamic ones: 

Cg(f) = Cf(f) + Cf"(f) 

The adiabatic part is defined as the response to the time-dependent source term g^^f{xY + w^^rn{x) in Eqs. (14), 
and is defined as the solution of the system of equations: 

° = ( - S^ + %)^f (^) + / (0%-?(^)Q'(^) + m^X? + Wj:^ra{x) (16) 



obtained by setting dCr/dt = in Eq. (14). It is formally solved as 



Cf{x) = j(fy jSqG^^{y- x)(g^^{yf + w^^^{y)) 



(17) 



where G^^{y — x) is a Green's function associated with Eq. (16). When introducing C^^ in Eqs. (13|l4), it appears 
only in integrals over momenta fig where the high-momentum contributions are the most significant. Tor sufficiently 
high momenta, we can use a simplified version of (17) 



Cf{x) 



5-g*(f)2 + w;g*^(f) 



% 



(18) 



This amounts to neglecting the effects of the trap, i.e. the integral term and the Laplace operator in Eq. (16). Using 
this expression, we finally obtain 



where 
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in—{x) 



ih 



dt 
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V{x) + g\^{x)\'^\^{x) + ^*{x){w^rn{x) + 

2V{x) + A - A' - i-i/i) *™(f ) + w'^{xf ^ 



(Pk 

(27r)3^^-fe 



^.c^i^ 



(Pk 

W) 



.wzciy-ix) 



(19) 
(20) 



w 



A' 



w 



f 



£k §1 



(27r)3 



{2nfEr 



(21) 
(22) 
(23) 



Unlike g and w , the coupling constants g and w are well-defined quantities. The constant g corresponds to the 
effective interaction strength between condensate atoms (one can show that g = g^ = Aiih^a/M where a is the s- 
wave scattering length associated to the potential U). The constant w is the effective coupling between condensate 
atoms and molecules, and one can show that w = w^. The quantity A' is a shift of the molecular energy due to the 
atomic-molecular coupling, commonly referred to as "light shift" in the case of an optical resonance. 
Equations (19), ([2Q|, and Eq. (17) along with the equation 



z^Cf "(f) = 



AM 



EuYtis) + J ^v,.ix)ciy-ix) - ^nl^cfix) 



(24) 



form a closed set of equations. Often, Cg deviates significantly from zero for small ^'s which lie in the Wigner's 
threshold law regime [27J. In this regime, one can make the approximation [31 J that g^;^ g and w^ ^ w. As a result, 
the whole set of equations is determined by the three molecular parameters g, w and A^ These quantities can be 
measured experimentally or calculated from the precise knowledge of the molecular potentials /7, Um and W. 



Solution of the equations for an instantaneous coupling 



further, we apply a local density approximation to C-. 



dyr 



In general, the resolution of Eqs. ( |19| , ( [2Q| ), and ( [24| ) is involved and requires numerical calculation. To proceed 

namely we neglect again the integral term and the laplacian 



in Eq. (24), 



^^§tCt{^) = E^,cly-{^ 



ih d 



{g-^{xf+W^^{x)^ 



(25) 



-ydyn 



SO that the spatial dependence of CJ^^ comes only from that of ^ and ^^. Thus, we make some error regarding the 
infiuence of the external trap V. This error should be small for times much smaller than the typical oscillation times 
in the trap. In the following, we will see that the interesting dynamics occurs at the microsecond timescale, while trap 
oscillations are typically around the millisecond timescale. This approximation is therefore justified. 



We are interested in the case where we start from a purely atomic condensate and suddenly turn on the coupling w 
to a molecular state. The initial conditions are 

*(f,0) = ^/M^ (26) 

^m{x,0) = 
Cf»(f,0) = 



where po{x) is an initial density profile. Integrating Eq. (25) gives 



This leads to 



Cf "(f , t) = I* ei^-^f-*) i- 1- (5*(f, r)2 + w<i>^{x, T)yT. 

I -^,Cf"{x) = l%{t-T)^[g^{x,Tf + w^Ux,T))dT (27) 

with the complex function 

The £ > is introduced to maintain convergence of the integral. It originates from the momentum dependence of i^^, 
which we ultimately neglected (which means that e is set to zero in the end). 



We now assume that the time integral in Eq. (27) is dominated by short-time contributions from ^^^, which 
leads to the Ansatz 



/ f{t-T)—(^g^l{x,Tf^w^^{x,T)^dT^af{t)w^,ri{x,t) 



(28) 



where a is a numerical factor to be determined. Inserting this result in Eqs. ( TopQ ), we finally get the closed set 
of equations ([3p| . 



To determine a, we consider the case where the transient shift and broadening is dominant in Eq. (|4]), which gives 



^m{x,t) 



wafit) 



Inserting this expression in Eq. (28), we find a = 7r/2. Interestingly, the same reasoning can be applied to systems 
of reduced dimensionality, which can be created with an external potential V{x) strongly confining the atoms in one 
or two directions |28l [29]. In this case, the motion of the atoms is frozen in the confined directions. This renormalises 
the coupling parameters g^ w and A' [30j, and the integration over momenta is reduced to 2 or 1 dimensions in all 
expressions. For a 2D-like system, we find 

f2D{t) = ^ge^^o^r(0,^^ot)^-^^log(^otM) for t^h/Eo 
(y-2D = 1 

where T is the generalised gamma function, and Eq is the zero-point energy due to the confinement. For a ID-like 
system, however, one finds aiD = 0, which indicates that the Ansatz is invalid in this case. 
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